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We use a simple analytic model to compute the angular correlation function of 
clusters identified in upcoming thermal SZ effect surveys. We then compute the ex- 
pected fraction of close pairs of clusters on the sky that are also close along the line 
of sight. We show how the expected number of cluster pairs as a function of red- 
shift is sensitive to the assumed biasing relation between the cluster and the mass 
distribution. We find that, in a ACDM model, the fraction of physically associated 
pairs is 70% for angular separations smaller than 20 arcmin and clusters with specific 
flux difference larger than 200 mJy at 143 GHz. The agreement of our analytic results 
with the Hubble volume TV-body simulations is satisfactory. These results quantify the 
feasibility of using SZ surveys to compile catalogues of superclusters at any redshifts. 

Key words: large scale structure of Universe - galaxies: clusters: general - cosmology: 
miscellaneous - methods: analytical 



1 INTRODUCTION 

Superclusters are regions with average mass overdensities larger than a few, on scales larger than a few megaparsecs; their 
cluster members show indications of intense dynamical activity and provide evidence that structures form hierarchically on 
supercluster scales (e.g. 



Bardelli, Zucca & Baldi 2001 



Rlnes et al. 2001 



Plionis fc Basilakos 2002) 



West 198£; 



Zucca et al 



Analyses of catalogues of superclusters in the local Universe (e.g. Bahcall & Soneira 1984 
199; ; Einasto et al. 2001), beyond the local supercluster ( |Oort 19*83 ), have focused on their shapes and their large scale 
distribution ( Kerscher 1998|). Sha pe statistics have been used for discriminating among the cosmological models (Basilakos 



Plionis & Rowan- Robinson 2001; Kolokotronis, Plionis & Basilakos 2002), whereas the total mass of superclusters can be 
used to estimate the density parameter fin from the mass-to-light ratio on large scales (Small et al. 1995). At high redshift, 



superclusters, when they exist (Postman et al. 2002), are difficult to identify, although preferred directions of radio emission 



from distant quasars and radio galaxies (West 1991) or unusual X-ray morphologies (Hashimoto et al. 2002) can provide 
evidence of their existence at z ~> 1. 



The peculiar velocities of clusters are enhanced in superclusters by non-linear effects (Zucca et al. 1993; Bahcall, Gramann 



fc Cen 1994). Colberg et al. (2000a) have indeed shown that the peculiar velocity of dark matter halos with massive neighbours 



in A-body simulations is ~ 40 — 50% larger than predicted by linear theory. In fact, linear theory predicts that the evolution 



of the peculiar velocities is independent of the local density, whereas Sheth & Diaferio (2001) have shown that, in iV-body 



simulations, evolved peculiar velocities of halos with the same mass are larger (smaller) in high (low) density regions. Moreover, 
the high velocity tail of the cluster peculiar velocity distribution can be a sensitive discriminator of cosmological models 
(Bahcall et al. 1994; Peel fc Knox 2002] ). In A-body simulations of representative volumes of the Universe, high density 



regions of a few megaparsec siz e contain two or more massive clusters; these superclusters can be easily spot in thermal 
Sunyaev fc Zel'dovich (1980; SZ|) effect surveys as nearby CMB te mperature decrements in the Ray le igh- Jeans limit and are 
the favorite regions for searches of enhanced kinematic SZ effect ( Diaferio, Sunyaev fc Nusser 2000 ). Confusion caused by 
primary Cosmic Microwave Background (CMB) anisotropies, thermal S Z and instrumental noise, compli cates the estimation 
of cluster peculiar velocities from measurements of the kinematic effect ( Aghanim, Gorski fc Puget 2001 ). The kinematic SZ 
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effect is difficult to detect directly; it is roughly an order of magnitude smaller than the thermal effect. One may proceed 
by looking for kinematic fluctuations in clusters selected according to thermal SZ measurements. However the most massive 
clusters that are responsible for the largest thermal fluctuations do not necessarily have the largest peculiar velocities. 

The thermal SZ effect, being independent of redshift, is the optimal tool for detecting clusters at high redshifts (Garlstrom 
2002|)] Several SZ surveys will produce catalogues of thousands clusters at any redshift in the next few years, e.g. SPT 
( |Carlstrom 2002[ ), AMI ( fKneissl et al. 200l[ ), AMiBA ( ]Lo et al. 2002| ), to mention a few. 

Clusters are not randomly distributed on the sky: their angular correlation function mirrors the clustering evolution of 
mass. The angular correlation function of SZ clusters was first computed by Cole & Kaiser (1988) who were interested in 
the non-Gaussianity of CMB temperature fluctuations caused by galaxy clusters. Komatsu & Kitayama (1999) show that the 
correlation in the cluster distribution enhances by ~ 20 — 30% the Poisson contribution to the CMB angular power spectrum 
at degree angular scales. Moscardini et al. (2002) compute the correlation function on the light-cone of clusters detectable 
by the Planck surveyor satellite and show that the correlation length of these SZ selected clusters is more sensitive to the 
physical properties of the ICM than to the cosmological parameters. 

In this paper we use a simple analytic model to compute the angular correlation function and the fraction of pairs that 
are close along the line of sight as well as on the sky. Our result shows that for viable cosmological models, this fraction is 
substantial. Therefore if we define superclusters as regions containing two or more close clusters, this result motivates the use 
of thermal SZ cluster surveys to identify superclusters at any redshift. We also compute the expected number count of cluster 
pairs within a given separation as a function of redshift. Unlike the number count of individual clusters, the pair number 
count is sensitive to the assumed form of the bias function between the cluster and the mass distribution. 

In section ^ we outline the basic equations; in section |^ we compute the expected number of clusters per solid angle in 
flux limited thermal SZ surveys, their angular correlation function and the number of physically associated pairs; we then 
apply our model to a full sky survey like the Planck surveyor. In section 3.5 we show how the fraction of physically associated 
pairs depend on the assumed biasing relation. We conclude in section W. 



2 BASIC EQUATIONS AND DEFINITIONS 

The interaction between the CMB radiation field and the hot cloud of electrons in the ICM plasma changes the specific 
intensity of the radiation. The difference in the specific intensity AI V at frequency v is (SZ) 

AI V = iog(x)y = i g(x)—^- / n e (r)T e (r)dl , (1) 
m e c^ J 

where Iq = 2(kT) 3 /(hpc) 2 , k is the Boltzmann constant, T — 2.725 K the present day CMB temperature ( Mather et al. 1999 ), 



hp the Planck constant, c the speed of light, <tt the Thomson cross section, m e the electron mass, x = hpv/kT, 



1 \ x ^ e * f e x + l \ 

and n e and T e the electron number density and temperature, respectively. The integral is over the line of sight I, r 2 = w 2 + l 2 ■ 
where w is the separation from the cluster center projected onto the sky, and we assume that the cluster is spherically 
symmetric. Here, we neglect the contribution of the kinematic effect. Let 

Jn e (r)T e (r)dV R 2 p f 

(Te> " = Jn e (r)dV ~NFJ "-WW* < 3 > 

be the electron density weighted mean temperature, where R p is the projected proper size of the cluster and N e ° the total 
number of electrons in the cluster. The variation of the CMB specific flux after passing through the cluster is 

Fu= / AIudQ » ALAn = i g(x) fc °" T 9 fg -Jrr^nM (4) 
J m e c 2 fi e m p d A A {z) 

where we have assumed an effective solid angle AQ = R 2 /d\(z), c?a(z) is the angular diameter distance, Nl ot — M f g / fj, e rn p , 
M is the total mass of the cluster, f g the gas fraction, /i e the mean molecular weight per electron, and m p the proton mass. 
Hereafter, we adopt f g = 0.06h~ 3 ^ 2 , which has been derived from a sample of 36 X-ray luminous clusters with redshift in 
the ranee 0.05 - 0.44 flEttori fc Fabian 1999[ ); here, for simplicity, we neglect (1) the observed variation, up to ~ 50%, of 



f g between clusters, and (2) the observed dependence of f g on redshift; note that this dependence is weaker in low-density 
Universes than in high-density Universes. Assuming a fully ionized ICM of hydrogen and helium with a helium mass fraction 
Y — 0.24, the mean molecular weight per electron is fi e = (1 — Y/2)~ x = 1.136. 

X-ray obse rvations indicate that the ICM temperature is roughly proportional to M 2//3 where M is the virial mass of 



the cluster (e.g. Finoguenov, Reiprich fc Bohringer 2001 Allen, Schmidt fc Fabian 2001), in agreement with arguments based 



on the virial theorem. Hydrodynamical simulations of cluster formation give similar slope for the T — M relation, but the 
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Figure 1. Upper panels: mass threshold corresponding to F™" 1 = 200, 50, 10 mjy (top to bottom) for v = 143 GHz (solid) and v = 353 
GHz (dashed lines). Lower panels: redshift cumulative cluster number counts corresponding to F™ ln = 10,50,200 mjy (top to bottom). 
Bold lines are the expected distribution (equation ph; thin lines are the actual number from the Hubble volume simulations. Dotted lines 
indicate the redshift limit of the lightcone outputs. 



temperatures are higher by a factor of two or so (e.g. Mathiesen fc Evrard 2001 ). Simulations might miss some important 
ingredient of the ICM physics; thus, we adopt aT-M relation in agreement with X-ray observations, i.e., 

2/3 



{T e ) n = Q.43[A e (z)H\z)/H$] 



2il/3 



M 



10 14 h~ 1 M Q 



keV 



(5) 



where h is the present day Hubble constant H in units of 100 kms" 1 Mpc -1 , H 2 (z) = Hi [fi (l+z) 3 +(l-flo-Ao)(l+2) 2 +Ao], 
and A c (z) is the mass density of the cluster, which has just collapsed at redshift z, in units of the critical density at that 
redshift (Bryan & Norman 1998). Note that here we assume valid the simple top-hat collapse model, although Voit (200C) 
shows that this assumption may be too simplistic when used to derive the evolution of the mass-temperature relation. We 
also assume that the X-ray temperature equals the electron density mean weighted temperature {T e ) n . 
The specific flux variation becomes 



F, = 9.33 x 10 4 (|)%(*)^f (^^\ 2 [A c (z)H 2 (z)/H^ 



M 



5/3 



mjy 



(6) 



SZ surveys will be able to detect only specific flux differences larger than a threshold F™ in . This constraint translates into a 
mass threshold Mth below which clusters remain undetected: 



M th (z,Fr n ) 



= 1.04 x 10" 



6/5 



- 3/5 M ( T^T^ ] [A c (z)H 2 (z)/H 2 ]-^ < / '"'"" 



(7) 



The top panels in figure ^| show Mth(z, F™ ln ) for two frequencies measurable by the High Frequency Instrument (HFI) 
of the Planck surveyor satellite in two cosmological models (see section |^ for details on the models). Except for an initial 
rapid increase, M t h varies slowly with redshift. This feature effectively makes a flux-limited SZ survey a mass-limited survey 



(Bartlett & Silk 1994) 
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3 STATISTICS OF THE CLUSTER DISTRIBUTION 



In this section, we compute the redshift distribution of individual clusters, the cluster angular correlation function w(9), and 
the redshift distribution of physically associated cluster pairs. 



We compare our analytic calculations with the lightcone outputs of the two Hubble volume simulations (Evrard et al 
2002 ; [Yoshida et al, 2001 ). These two simulations correspond to a ACDM and a rCDM cosmology with (fio, Ao, h, as) = 



(0.3, 0.7, 0.7.0.9) and (1, 0, 0.5, 0.6), respectively. For both models, the power spectrum shape parameter is F = 0.21 (Frenk et 
al. 20(Q; see Evrard et al. 2002 for further details) . Using the lightcone outputs has the advantage of simulating a more realistic 
observation of a patch of the sky when compared to the usual technique of stacking simulation boxes at different redshifts 
(e.g. |5caramclla, Cen fe Ostriker 1993| ; |da Silva et al. 2000| ; |5eljak, Burwell fc Pen 2000|; |Springel, White & Hernquist 2001 



Zhang, Pen fc Wang 2002j ). Dark matter halos in the A^-body simulations have been identified with a spherical overdensity 



technique. We assume that these halos do correspond to galaxy clusters in the appropriate mass range, although comparisons 
with the APM cluster catalogue show slightly different clustering properties (Colberg et al. 2000b). 



3.1 Number count of individual clusters 



The number count rj per unit solid angle of clusters as microwave sources, as compared to X-ray sources, was first computed 
by Korolev, Sunyaev & Yakubsev (1986); many authors, which were interested in different aspects of r/, have computed it 
ever since (see e.g. Bartlett 2001 and references therein). The number count rj enters the calculation of the angular correlation 
function (section 3.2). It is sensitive to the properties of the ICM (Holder & Carlstrom 2001; Benson et al. 2002]) and to the 



2001 



cosmological parameters (Barbosa et al. 1996; Kneissl et al. 2001; Haiman, Mohr & Holder 2001; Holder, Haiman & Mohr 



Benson, Reichardt & Kamionkowski 2002). The number count also depends on the parameters of the survey, namely 



its sensitivity (or its flux limit FJ m ) and its angular resolution. The cumulative number count of clusters, at redshift greater 
than z, detected per unit solid angle is 



r\(z, F™ ln )dSl = dfi 



rf 2 (zi) 
H(zi) 



dzi 



dn(M,zi] 
dM 



dM 



(8) 



'AI th ( zl ,F»'") 

where d(z) — (1 + z)dA(z) is the distance measure, and dn(M,z)/dM is the comoving number density of clusters with mass in 
the range (M, M + dM) at redshift z. We have assumed that the detector is able to measure the total specific flux difference 
F v of each cluster. When the angular resolution of the survey is included in the analysis, the cluster number count can drop 
significantly (Bartlett 200C). 



In the Hubble volume simulations, the dark matter halo mass function is well fitted by the formula (Jenkins et al. 2001) 
dnj(M,z)_ dM _A(p) 



M 



dM 

where {A, B, e) — 
mean density, a e g 

In a' 1 (M) = - hi (lis +alnM + 6(lnM) 



a e ff(M) exp[-| In a~ 1 (M) + B\ e ]dM 



O) 



(0.22, 0.73, 3.86) and (0.27, 0.65, 3.77) for the ACDM and tCDM model, respectively, (p) is the comoving 
= din a- 1 {M) j din M, 

(10) 



and (cri5, a, b) = (0.578,0.281,0.0123) and (0.527,0.267,0.0122) for the ACDM and rCDM model, respectively. 

The lower panels of figure [l] show the redshift cumulative cluster number count for the two models (thin lines) compared 
with the expectation of equation (^|). We use the NO lightcone outputs of the Hubble volume simulations, which cover -it/ 2 sr 
and have depth z = 1.46 (ACDM) and 1.25 (rCDM), as indicated by the dotted lines. The disagreement at low J 71 ™ 111 and at 
high redshift is mainly due to the mass resolution limit. 



3.2 The angular correlation function 



To compute w(8), we need to know the three dimensional two-point correlation function £(ri2, Mi, M2, Zi, Z2) for clusters of 
mass Mi, M2 at redshift Z\, z% and comoving separation ri2. When z\ — 22 = z, we can use the simple biasing scheme 



Z(ri2,Mi,M 2 ,z) = &st(Mi,z)&st(M2,z)£ 



where flSheth fc Tormen 1999| ) 

— 2 1 . 2p 



(11) 



(12) 



b S T(M,2;) = 1 + . . + T . . — jr-y, 
S c (z) S c (z)[l + (av 2 )P\ 

a = 0.707, p = 0.3, v = S c (z)/a(M, z), <j{M, z) is the rms of the mass density field smoothed on scale R = (3M/47r(p)) 1 ^ 3 , (p) 
is the mean comoving mass density, and 5 c (z) is the li near overdensity extrapo lated at redshift z for a spherical perturbation 
which has just collapsed (see Kitayama fc Suto 1996 or Nakamura fc Suto 1997 for handy fitting formulae to these quantities); 
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£ m (r, 2) is the correlation function of dark matter at redshift z. Equation ( |l^ ) is strictly valid on large scales. Because we are 
interested in close pairs of clusters, we need to consider a scale dependent bias factor b(r, M, z) (e. g. Yoshikawa et al. 2001). 
For clusters at different redshifts Z\ and Z2, we can modify equation ( |ll| ) to 

£(ri2 7 M 1 ,M 2 ,zi,Z2) = b(ri2,M 1 ,zi)b(ri2,M 2 ,Z2)Zm(ri2,zi) (13) 

where we break the symmetry in z\ and 22 by "sitting" on the cluster at z\. Equation ( |l3| ) is consistent with equation (|l"lj), 
because clusters at different epochs Zi, Z2 are not correlated and €,{ri 2 ) goes to zero with increasing ri2 quickly enough that 
£(^12) 7^ only when z\ ~ 22. Alternative choices of the redshift at which £ m can be computed are discussed in 
(1997|. 



Porciani 



14.2.7; Osmer 1981) 



The comoving relative separation rn in a space described by the Robertson- Walker metric is (Weinberg 1972, equation 

(14) 



rf 2 (zi, Z2, 8) = d 2 {z 1 )G 2 (z 1 ,z 2 , 9) + d 2 {z 2 ) - 2d(zi)d(z 2 )G(zi, z 2 , 8) cos0 



where 8 is the angular separation of the two objects on the sky, 

d ( 22 ) „„„n 



G(zi,Z2,6) = y / l-Kd 2 {z 2 )+ (l - Vl-^ 2 («i)) 



d(zi) 



(15) 



and re = (Ho / c) 2 (S7q + Ao — 1). This equation is not symmetric in zi and 22 for a non flat Universe (k / 0). However, if 
we deal with distances smaller than the curvature radius of the Universe l/^/fre|, we have 1/re 3> d 2 (z\) and 3> d 2 (z 2 ), and 
G ~ 1; this is the case of interest, because we expect £(^12) ~ at large separation. When z\ — z 2 — z, r\ 2 ~ 2d 2 (l — cos8) 
for K,d 2 ~ 0, and we correctly recover 8 = ri 2 /d(z) for small 8. 



et al. 



1997 



To compute the angular correlation function, we follow the classical argument that leads to Limber's equation (Matarrese 

nding two clusters with diffei 

dM 2 [1 + £(r 12 , Mi , M 2 , zi , z 2 



Moscardini et al. 2000). Consider the probability of finding two clusters with different mass at different redshift 

dn(M 2 ,z 2 ) 



C(r 1 2,M 1 ,M2,zi,z 2 )dAidA2 = dVidV 2 dn \,i' Zl ' dMi 



(16) 



dMi dM 2 

where dA t = dMidzidfli, dfli is the infinitesimal solid angle and dVi = cd 2 (zi)dVL i dz i j'H(zi) is the infinitesimal comoving 
volume. The probability of detecting two clusters at angular separation 8 is proportional to 

ip(8)dQ, 1 dQ, 2 = df2idf2 2 / C(»"i2, Mi, M 2 , zi, z 2 )dM 1 dM 2 dz 1 dz 2 

(17) 



= d^id^2Vo{Ff U ) [l + Fr n )] 
where ri (F^) = 77(0, *T n ), and 
Ho(e,Fr n ) 



w(8,Fr n 



(18) 



is the angular correlation function we are seeking. All the clustering and biasing information is in the function 

roc rp-(y \ 

• — 1/ r\ 7-,min\ 2 / / j-,min\ " 1*1 } 7 

~{z,8,F„ ) = c n(zi,F v ) ; ' dzi x 

d 2 (2 2 ) 



■ f n(z 2 , -F™ m )& c ff(ri2, Zi, F™ ln )b oH (r 12 , z 2 , F™ ln )£ m {r 12 , 21) 

I H{Z2) 



dz 2 



where 

n[z,F v ) = 



dn(M, z) 



M th ( Z ,F»'°) 
1 1 



dM 



dM 



b eS (r,z,Fr n )=n- L (z,Fri 



and E (8,Fr n ) = E(0,8,Fr" 



,. .dn(M,z) r 

6(r, M, z) — \-r7— L dM 

K ' ' ; dM 



(19) 

(20) 
(21) 



We use the fitting formula given by Peacock & Dodds (1996) to compute the correlation function of dark matter f m (n2, 2) 
in the non-linear regime. For the effective scale dependent bias factor b e « we adopt the following fitting formula. 



b e a(r,z, F™ in ) = &sT, c ff(z, F™ m )[l + bsT, c ff(2, F™' n )a(r, 2)] ' 35 , 
where er(r, 2) is the mass variance smoothed over the top-hat radius r, and 
Jr th (*,F--) 6 st(M, z)[drv(M, z)/dM]dM 



&ST,eff(z, -F 1 ™) = 



[dnj(M,2)/dM]dM 



(22) 



(23) 
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Figure 2. Angular correlation function of SZ clusters with F™ m > 200, 50, 10 mjy (top to bottom ) for v = 143 GHz (left panels) and 
v = 353 GHz (right panels). Bold lines are the expected angular correlation function (equation |l8|); filled squares are the correlation 



function from the Hubble volume lightcones. 



Figure ^| shows that this choice of the biasing relation yields an excellent agreement between the model and the angular 
correlation function derived from the pair count of SZ clusters in the iV-body simulations. The functional form of the biasing 
relation (equation E2J) was first suggested by Hamana et al. (20011) , with a power 0.15 rather than 0.35. Hamana et al's relation 
was only checked against the redshift zero output. When we compare the analytic modelling with the lightcone output, we 
include a broad redshift range, and the agreement shown in figure |^ indicates that a larger power seems more appropriate. 



3.3 The number of physically associated cluster pairs 

The fraction of pairs with angular separation smaller than 9 whose cluster members have comoving separation along the line 
of sight smaller than r max is proportional to an expression similar to equation ( p^[ ) where the integral over 22 is limited to 
the interval set by r max . Because £(^12) ~ at large ri2, when r max is of order the cluster correlation length, we can estimate 
this fraction as 



¥*>(< 0) = 



/„' w{9 1 )9 1 d0 1 
/ o e [l + u>(0i)]0id0i 



(24) 



Thus, ifio(< ff) also indicates the fraction of close pairs of clusters which are roughly at the same redshift. 

When comparing <fio(< ff) with the actual pair count in the TV-body simulations, we need to quantify what we mean 
by close along the line of sight. We compute ipo{< ff) with w(0) derived from the simulations and we choose the maximum 
comoving pair separation that provides a reasonable agreement between </'o(< 0) and the actual count of close pairs. The 
comoving volume of the rCDM model is (2/3) 3 smaller than the comoving volume of the ACDM model, but the two volumes 
contain the same total mass; we thus choose, for the two models, two different maximum comoving separations whose ratio 
is 2/3: therefore, we have the same mass within the volumes defined by these separations. Figure ^ shows that the agreement 
between <po(< 9) derived from the simulations (bold lines) and the actual pair count (squares) with cluster member three- 
dimensional comoving separation smaller than 16 (10.7) ft -1 Mpc for the ACDM (rCDM) model is satisfactory. Dotted lines 
show <po(< 9) computed with the analytic model. The agreement between the simulations and the analytic model is excellent. 

Figure |3| shows that the analytic model predicts that more than ~ 70% of the pairs with angular separation 9 < 20 
arcmin are physical associations in the shallow survey F™ ln > 200 mjy at v — 143 GHz. Because of the larger number of 
clusters per square degree, the fraction of pairs close in three dimensions drops to ~ 30% in the deeper survey F™ 111 > 50 
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Figure 3. Fraction of cluster pairs with comoving separation < 16.0 (10.7) h^ 1 Mpc for the ACDM (rCDM) model (squares). Upper 
(lower) curves are for F™ ln > 200 (50) mjy. Left (right) panels are for v = 143 (353) GHz. Bold lines correspond to <gnf< 8) where 
w{0) is computed directly from the TV-body simulations. Dotted lines are the corresponding analytic functions (equation K4J) . Error bars 
assume Poisson statistics. 



mjy; this fraction is still substantial. Moscardini et al. (2002j ) show that, because the properties of the ICM affect the cluster 



mass threshold, for a fixed specific flux threshold F™ ln , the correlation length ro measured for SZ clusters is more sensitive 
to the ICM properties than to the cosmological parameters. According to their results, however, our estimated fractions of 
physically associated pairs can change by 30% at most. 

The cumulative redshift distribution of the fraction of physically associated cluster pairs with angular separation in the 
range (9i, 62) is 

<p{z,0u0a) = — g . (25) 

Je' lVo(Fr n ) + S o (0, F^)]9d9 

Figure ^| compares the model with the simulations. The agreement is satisfactory. Shallower surveys will detect only massive 
clusters which form at late times. Therefore most pairs of close clusters are at moderately low redshift. 

To compute the number of pairs with separation 9 on the sky, consider the probability rjg[l + w(9)]dQidfl 2 of having 
such a pair, where dSli is the infinitesimal solid angle around each pair member and the angle between the versors pointing 
to dili and dtt2 is 9. Without loss of generality, we can fix the versor of dfii coincident with the 2-axis and write dQ.2 = 
d</>2dcos#2 = d(j>2dcos8. Integrating over d£lid(f>2 yields J dQ\d£l2 = 87rsinf)df9. All over the sky, the number of pairs with 
separation in the range (6, 9 + d9) is thus h p (8)d9 = 87r 2 r;o[l + w(6)] sin9d9/(2 + (w)), where (w) = J* sin9w(9)d9 and the 
relation is normalized such that it yields a total of (Anr/o) 2 /2 pairs at any separation. In a survey with area Af2, the number 
of pairs with angular separation in the range {9\, 62) is trivially (Afi/47r) J® 2 h p (9)d9. We use this relation in the next section. 



3.4 A survey example 

To get a better idea of the sample sizes, consider a full sky survey, like the one that will be performed, for example, by the 
Planck surveyor, with a specific flux difference limit 100 mjy at 143 GHz. In a ACDM (rCDM) model, this survey would 
detected ~ 5500 (17400) clusters and ~ 155 (1150) pairs with angular separation in the range 8-16 arcmin. Of these pairs, 
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Figure 4. Cumulative redshift distribution of the fraction of cluster pairs with comoving separation < 16.0 (10.7) h~ 1 Mpc, F™' n > 50 
mjy, and angular separation in the range (8, 16), (16, 32) and (32, 64) arcmin (top to bottom) for the ACDM (tCDM) TV-body models 
(thin lines). Left (right) panels are for v = 143 (353) GHz. Bold lines show the corresponding analytic distributions <p(z, 6\, 82) (equation 



83% (53%) are also close along the line of sight. It means that measuring the redshift of 155 x 2 (1150 x 2) clusters in close 
pairs would yield 128 (614) physically associated pairs and 54 (1052) "isolated" clusters. These numbers are computed with 
the analytic model described above. The differential redshift distributions n c (z) of the 310 (2300) clusters and n p (z) of the 
128 (614) pairs are shown in figure |B|. 

Extracting the cosmology information from the cluster redshift distribution is not trivial. The redshift distribution of 
clusters depends on the cosmological parameters in a complicated way, through the geometry of the Universe and the cluster 
mass function (equation ^) . The redshift distribution of pairs introduces further complications due to the cluster bias function 
and the non linear evolution of the dark matter correlation function (equation [l9| ). So an estimation of the cosmological 
parameters with a traditional least squares fitting procedure is very cumbersome. A simpler procedure is to compare the 
observed distributions with those predicted by a few fiducial models. Therefore, we perform the following Monte Carlo 
simulation. First, by assuming a cosmological model, we use our analytical model to compute the expected number of clusters 
n c (z) and pairs of clusters n p (z) in a given redshift bin. We then generate samples of clusters and cluster pairs from Poisson 
distributions with mean values n c (z) and n P (z), respectively. 

Finally we compute 

zbin 

where the sum is over the redshift bins, the indices i and j refer to the cosmological model (e.g., ACDM or rCDM), rij is the 
variable drawn from model j and is the expected mean value in model i. 

The lower panels of figure |^ show the distribution of the \ 2 f° r 10 4 Monte Carlo realizations drawn from the redshift 
distributions shown in the upper panels. Each random sample (either of clusters or pairs of clusters) drawn from a model j is 
a point on the plane x\\j ~~ Xan m figure |^. Both the cluster and pair samples would give significant large \ 2 when compared 
to the wrong cosmological model. 



3.5 Sensitivity on the biasing relation 

The predicted number of physically associated cluster pairs depends on the assumed biasing relation between the cluster and 
the mass distribution. In the previous sections, we have assumed a time and scale dependent bias factor (equation To 
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Figure 5. Upper panels: differential redshift distribution n p (z) of physically associated pairs of clusters with angular separation in the 
range 8—16 arcmin for a full sky survey with F™ m = 100 mjy and v = 143 GHz (solid histograms); differential redshift distribution 
n c (z) of all the clusters in pairs with angular separation in the range 8—16 arcmin, irrespective of their physical association, in the same 
survey (dashed histograms). Lower panels: probability density of the reduced x 2 of Monte Carlo realizations drawn from the tCDM 
model (left panel) and the ACDM model (right panel) computed with respect to each model. Solid (dashed) contours are for the pair 
(cluster) samples. Contour levels are (0.61,0.13,0.01) times the maximum probability density. They roughly correspond to the 68%, 95% 
and 99% confidence levels. 



investigate how sensitive the fraction of pairs is on this assumption, we consider two extreme models: a scale independent and 
a time independent bias factors. Namely, 

b no - r (z, F™ n ) = bsT.cff (2, iC") (27) 

and 

b no _ z (r,Fr n ) = fe ST ,cff(0,Fr n )[l + fe ST , cff (0,Fr>(r,0)] ' 35 . (28) 

Figure ^| shows our fiducial biasing relation b c s (r, z, _F™ m ); in this figure, the bottom lines of the ACDM and rCDM 
models show the time independent bias factor b no - z {r, F™ in ). Figure § show the scale independent bias factor b no - r (z, _F™ m ). 

Figure ^ shows that neglecting the scale dependence considerably suppresses the angular correlation function (dotted 
lines). Because clusters form at relatively low redshift, neglecting the time dependence is less dramatic (dashed lines): in fact, 
the number of physically associated pairs is substantially underestimated only at z > 0.2 (bottom panels). 



To quantify how significant the difference between these biasing relations is, we proceed as in section 3.4. We draw 10 4 
Monte Carlo realizations from the redshift distributions of the physically associated cluster pairs shown in the bottom right 
panel of figure |^. We compute the variable x|u (equation E^), where now i and j stay for our fiducial (equation |2^), scale 
independent (equation |27|) or time independent (equation |28|) biasing relations. 

Our results are shown in figure |^. The difference between our fiducial model and the scale independent model generates 
large x 2 values when compared to each other (top panels) . The difference between our fiducial model and the time independent 
model, although substantially smaller, is still well measurable (bottom panels). 



4 DISCUSSION AND CONCLUSION 

One of the main purposes of SZ cluster surveys is to constrain the cosmological model with the cluster redshift distribution. 
In fact, the redshift distribution of SZ selected clusters depends on the density parameter f2o, the mass power spectrum, and 
the physical properties of the ICM. Therefore one hopes that catalogues of SZ clusters can constrain the cosmological models 
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Figure 6. Effective bias (equation |22j ) of SZ selected clusters in a survey with F™ m = 50 mjy and v = 143 GHz. Solid (dashed) lines 
arc for the ACDM (rCDM) model. Curves are for z = 0.01,0.02,0.06,0.14,0.34,0.82,2.0,4.8 (bottom to top). 
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Figure 7. Effective scale independent bias (equation p7[ ) of SZ selected clusters in a survey with _F,™ m = 50 mjy and u = 143 GHz. 
Solid (dashed) line is for the ACDM (tCDM) model. 
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Figure 8. Sensitivity of our results on the biasing relation in the ACDM model for a survey with _F™ m = 50 mjy and v = 143 GHz. 
Solid, dashed and dotted lines refer to our fiducial (equation p2| ), time independent (equation and scale independent (equation |2t| ) 
biasing relations, respectively. The bottom right panel shows the redshift distribution of physically associated cluster pairs with angular 
separation in the range 8-16 arcmin in a full sky survey. 
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Figure 9. Values of the x 2 corresponding to the bottom right panel of figure |j. See text for more information. Contour levels are 
(0.61,0.13,0.01) times the maximum probability density. They roughly correspond to the 68%, 95% and 99% confidence levels. 
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despite the degeneracies between some parameters such as the normalization of the power spectrum as and flo (Barbosa et 
19§; 



Kneissl et al. 2001 



Holder, Haiman & Mohr 2001; Benson, Reichardt & Kamionkowski 2002). 



al 

However, for the very same reason why SZ cluster surveys are so attractive, their catalogues will not have any accurate 
information on the cluster distances, although the morphology of the cluster can provide a rough estimate of its redshift (Diego 
et al. 2002). The redshift distribution of clusters can be determined with follow up optical observations, CO line measurement 
in galaxy cluster members, or high resolution X-ray spectroscopy (e.g. Bleeker & Mendez 2002). The number of objects in 
SZ cluster surveys depends on the sky coverage, instrumental sensitivity, and observational bands. In an all sky surveys, this 
number can be of order several thousands. Of course, only the redshifts of a subsample of these clusters is sufficient to provide 
a robust estimate of the redshift distribution. How should one choose this cluster subsample? We suggest to measure the 
redshift of clusters in close pairs. Because the probability that clusters close on the sky are also close along the line of sight is 
substantial, this strategy will provide, at the same time, the redshift distribution of clusters and a catalogue of superclusters. 

To show that the compilation of a catalogue of superclusters from SZ cluster surveys is indeed feasible, we have computed 
the angular correlation function of SZ selected clusters for two cosmological models. We have then computed how many 
cluster pairs, close on the sky, are also close along the line of sight. For example, for clusters with specific flux difference 
larger than 200 mjy at 143 GHz in a ACDM model, 70% of the pairs with angular separations smaller than 20 arcmin are 
physical associations. Because an SZ supercluster catalogue will contain superclusters at any redshifts thanks to the redshift 
independence of the SZ effect, the catalogue will be of great relevance for follow up studies of the evolution of galaxy and 
cosmic structures in high density regions. Moreover, from the redshift distribution of physically associated cluster pairs, we 
can extract information on the biasing relation between the cluster and the mass distribution. 

To compute the number of detectable clusters, we have used the point source approximation, where the survey can detect 
the total cluster radio flux, but it cannot resolve the cluster structure. This approximation is thus independent of the cluster 
profile. It is realistic for poor resolution surveys (for example the Planck surveyor will perform a survey with a beam FWHM 
of 8 arcmin at 143 GHz), but of course unrealistic for high resolution bolometric or interferometric ground-based surveys. 
However, resolving the clusters has the net effect of increasing the mass threshold and thus decreasing the number of detected 
sources (Bartlett 200C). In other words, resolving the cluster structure is similar to considering our results for shallower 
surveys: the cluster minimum mass is larger and the fraction of physically associated pairs increases. 
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